Non-local correlation effects and metal-insulator transition in the s-d 

exchange model 



J. Sweep, ^ A. N. Rubtsov,^ and M. I. Katsnelson^ 

^ Radboud University Ntjmegen, Institute for Molecules and Materials, 
Heyendaalseweg 135, NL-6525AJ Nijmegen, The Netherlands 
^Department of Physics, Moscow State University, 119992 Moscow, Russia 

(Dated: April 1, 2013) 

The metal-insulator transition within the s-d exchange model is studied by the Dynamical Mean 
Field Theory and Dual Fermion approaches. The latter takes into account nonlocal correlation effects 
which are shown to be essential. In particular, the critical values of the s-d exchange coupling 
constant for these two methods turn out to be different by a factor of more than two (for the 
case of square lattice and spin 1/2 localized electrons). The calculations were performed using a 
Continuous- Time Quantum Monte Carlo method. The difference of the quantum spin- 1/2 case and 
the classical-spin case is discussed. For the quantum case the sign of s-d exchange coupling constant 
is relevant which demonstrates an importance of the Kondo effect. 
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I. INTRODUCTION 

The motivation to study the metal-insulator 
transition in the s-d exchange, or Vonsovsky- 
Zener, model [U E] (now it is frequently called 
Kondo lattice model) is given by previous successes 
obtained studying the metal-insulator transition 
(MIT) in the Hubbard model [3 . The s-d exchange 
model is complementary to the Hubbard model. 
While in the Hubbard model the interacting elec- 
trons responsible for magnetism and for electron 
conductivity are located in the same band, in the 
s-d exchange model they are located in different 
bands. In one band the electrons are delocalized or 
itinerant, therefore these electrons are called con- 
duction electrons. The other band contains elec- 
trons that are localized on lattice sites. The type 
of the interaction is also different between the two 
models. The Hubbard model contains an on-site 
Coulomb interaction modeled by a single energy 
parameter U. The s-d exchange model contains 
a spin-spin interaction between the localized and 
conduction electrons with a coupling constant (s-d 
exchange parameter) J. The s-d exchange model 
is used in the description of rare earths, in which 
4f-electrons have a very localized nature [4l. The 
description of colossal magnetoresistance materi- 
als such as (La,Ca,Sr)MnO compounds sometimes 
is based on the s-d exchange model assuming that 
the t2g electrons are rather localized and the Cg 
electrons are itinerant |S]. 

Whereas for relatively weakly correlated sys- 
tems the crystal with half-filled energy bands is 
a metal, for strong enough correlation the sys- 
tem should be an insulator, with each electron sit- 
ting in a localized state associated with a given 
site. This corresponds to the splitting of the en- 
ergy bands into Hubbard subbands separated (for 



a large enough interaction parameter U) by a gap 
[BJ |7] • In the s-d exchange model the parameter U 
is replaced by J • S* where S is the spin of localized 
electrons. The Hubbard-III approximation [8, 9] 
provides a qualitative description of the metal- 
insulator transition for both models. This descrip- 
tion however does not take into account all quan- 
tum effects; for the case of s-d exchange model, it 
turns out to correspond to the classical spin limit 
-> oo jni- The Dynamical Mean Field Theory 
(dmft) [lIT was a real breakthrough in this respect 
dramatically improving the description of metal- 
lic phase near the transition. Whereas within the 
Hubbard-III approximation the main effect of ap- 
proaching the MIT is the growth of damping of 
electron states at the Fermi energy, together with 
appearance of shoulders as the precursors of future 
Hubbard bands (two-peak structure), the dmft 
gives the three-peak structure, with the "Kondo 
peak" near the Fermi energy, that is, the growth of 
effective mass instead of the growth of the damp- 
ing. The origin of this peak is the formation of 
a resonance due to electron scattering by quan- 
tum spins |10| . The dynamical mean-field theory is 
formally exact in the limit of infinite space dimen- 
sionality whereas the Hubbard-III approximation 
is formally exact (within the s-d exchange model) 
in the limit of infinite space dimensionality and in- 
finite, that is, classical, spin [9 . In particular, in 
this limit there is no difference between the case 
of ferromagnetic (J < 0) and antiferromagnetic 
(J > 0) s-d exchange coupling; the Kondo reso- 
nance exists only in the latter case (the detailed 
discussion of the Kondo effect and related physics 
can be found in Ref iTT)l . 

The DMFT is the best local approximation, that 
is, it corresponds to an optimal self-consistent 
choice of the momentum-independent electron self- 



energy. For real two- and three-dimensional sys- 
tems it is just an approximation, which is not quite 
controllable. The Dual Fermi on formalism |12| - 
114) allows to construct a diagrammatic expansion 
for the non-local correlation effects starting from 
the DMFT Green's function as a zeroth-order ap- 
proximation. The differences between the Dual 
Fermion results and the dmft results were inves- 
tigated recently [TS] for the Hubbard model on tri- 
angular lattice, and turned out to be very essen- 
tial. The present paper discusses the differences 
between Dual Fermions and dmft for the metal- 
insulator transition in the s-d exchange model. 

The MIT in the s-d exchange model in the clas- 
sical limit 5 — > oo has been studied in Ref 9 In 
this work we consider a quantum s-d exchange 
model at a two-dimensional (2D) square lattice, 
for 5* = I where we expect quantum effects to be 
maximal. We study the role of non-local correla- 
tion effects by comparison of the results of dmft 
and Dual Fermion calculations. We show that the 
DMFT strongly overestimates the critical value of 



the coupling constant for the metal-insulator tran- 
sition Jc- We also show that in the s-d exchange 
model quantum effects are indeed relevant, since 
the results for positive and negative values of the 
coupling constant J are essentially different. 



II. HAMILTONIAN 

The s-d exchange model contains electrons that 
are localized on lattice sites, and conduction elec- 
trons that can hop between different sites 
Let us introduce their annihilation operators Cia-, 
a —'[, 4, is the spin projection. The annihilation 
operators for the localized electrons will be desig- 
nated dia-. The Hamiltonian contains (in order of 
appearance in eq. ([ij) a hopping term, a Hubbard 
term that is included for calculational purposes, 
the s-d exchange interaction, and the chemical po- 
tential fi term: 
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N is the number of conduction electrons and nu„ — d\^dia. The spin-spin interaction term reads: 



The Hubbard interaction is included to (computa- 
tionally) ensure the half-filling of the lattice sites 
(one localized electron per site) . Strictly speaking, 
it is reached only in the limit U ^ 00. We will do 
the calculations for a large but finite U. With this 
Hubbard-term and U ^ J and f3U ^ 1 half-filling 
is achieved to sufficient extent. 



III. METHOD 

Our goal was to study the density of conduction- 
electron states (DoS) of the s-d exchange model at 
different coupling constants J. The calculations 
were carried out at the Matsubara axis (imaginary 
frequencies) and, to obtain the densities of states 
we performed an analytical continuation of the 
Green's functions on the real axis, using the Max- 
imum Entropy method |16l[T7] . The Green's func- 
tions we obtained self-consistently in dmft versus 
the Dual Fermion approach. In short, both dmft 
and Dual Fermion utilize an effective-medium con- 
cept, so that they deal with the impurity model 



described by the action S^^p = Sat + A^cj^^c^j-, 
where Sat is the action of an isolated lattice site, 
and A is a hybridization by the effective bath. 
Frequency-dependence of the hybridization allows 
to handle the effects of local fluctuations, that un- 
derlie the Kondo physics. On the other hand, this 
w-dependence makes the impurity problem non- 
hamiltonian and therefore requires a usage of nu- 
merical solvers to obtain its Green function g^^ and 
two-particle vertex 7^^^ (the later is required for 
Dual Fermions). As a numerical solver, we used 
the interaction-expansion continuous-time quan- 
tum Monte Carlo procedure [HI [T^. The Dual 
fermion lattice Green's function equals (spin and 
orbital indices are omitted) 



where is a dispersion law arising from the hop- 
ping tik, and E is the self-energy function of the 
dual ensemble obtained in certain diagrammatic 
approximation. Vertices of dual diagrams are ver- 
tex parts of the impurity problem. DMFT can 
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be represented as the zeroth-order Dual Fermion 
scheme with S = 0. In this approximation, the 
lines (propagators) in dual diagrams have a simple 
meaning, as they are equal to the non-local part of 
the Green's function, Gj^'jf^'^ = G^,fe— Gtj,r=o (for 

a general case, G and G are also exactly related, al- 
though the relationship is more complicated). One 
can easily see that DMFT describes purely local 
contribution to the lattice self-energy (and is actu- 
ally shown to correspond to the summation of all 
local diagrams), whereas Dual Fermion diagrams 
describe non-local effects. More technical detail 
can be found in Refs. |12l [TB] . 

The hybridization function was obtained 
self-consistently from the condition that the local 
part of the dual propagator vanishes, G^^r=o = 0. 
This eliminates the first-order Hartree dual dia- 
gram, and thus the leading correction to E is given 
by the second-order (see diagram (b) and formula 
(20) of Ref.[13]), 



1 

2/32 



E 



(4) n n n 



(4) 

(summation over spins is also implied). 

The last expression is an approximation used for 
the Dual Fermion calculations in the present pa- 
per. The calculations were performed for the para- 
magnetic state of the Hamiltonian ([T|) on a square 
lattice with the nearest-neighbor hopping parame- 
ter t = 0.25 (which means a half-bandwidth equal 
to 1), we have chosen U — 5 and the inverse tem- 
perature /3 — 10. Dual fermion diagrams were cal- 
culated at a 24 X 24 slab in momentum space. The 
final accurate calculations took about two weeks of 
CPU load per J-value in the dmft case and about 
10% more in the Dual Fermion case. The amount 
of self-consistency loops was roughly 10-20 while 
the number of Monte Carlo steps was roughly dou- 
bled with each loop in the last 6 loops until it 
reached 8 • lO^. 



IV. RESULTS 



Critical coupling constant for 
metal-insulator transition 



The MIT for the s-d exchange model is visible 
in Fig. [T] We have quoted values of the coupling 
constant just around the transition. Although the 
temperature is quite high and the Hubbard- [/ = 5 
we still see a qualitatively correct picture of the 
densities of states, accurate enough to estimate a 
reasonable value of the coupling constant J for the 
metal-insulator transition. We consider the state 
to be conducting if there is a local maximum at the 



Fermi level and insulating if there is a local mini- 
mum at the Fermi level, because the temperature 
is so high. 

The choice of L/ = 5 in the Quantum Monte 
Carlo calculations is justified because in our cal- 
culations (3U = 50 1 and U/J > 5 so that 
half-filling is achieved to sufficient extent. 

In Table |T] we compare several different calcula- 
tional methods that were used to obtain a criti- 
cal value of the coupling constant for the metal- 
insulator transition. 



B. Difference for positive and negative 
coupling constant 

The difference in behavior of the DoS between 
positive and negative values of the coupling con- 
stant is displayed in Figures [2] and |3] 

In Figure |2] we display the results for \J\ = 0.4 
calculated within both dmft and Dual Fermions. 
^,^The different colors represent positive and nega- 
tive J- values. Stars correspond as before to the 
DMFT-calculation while the Dual Fermion result 
is drawn with lines. In this picture we look at 
the general shape of the curves. We see that they 
are not the same for positive and negative J and 
that incorporating non-local corrections through 
Dual Fermions changes the results: the difference 
is more pronounced in the Dual Fermion case. 

In Figure |3] we display the results for \J\ = 
0.7, also calculated using dmft (stars) and Dual 
Fermions (lines). The different colors again repre- 
sent the positive and negative J- values. The gen- 
eral shape of the curves is much more different 
between positive and negative J for \J\ = 0.7 than 
for I J| = 0.4. 



V. CONCLUSIONS 

In this work we justify two claims. First, the 
Dual Fermion formalism description of the metal- 
insulator transition in the s-d exchange model 
to date predicts a critical value of the coupling 
constant Jc « 0.3. This value differs from the 
DMFT-result by a factor of more than 2 and corre- 
sponds nicely to the result obtained in the other 
work using a non-local correlations approach [20J. 
The different results between various calculational 
schemes are summarized in table |T1 

Second, there is a clear qualitative difference be- 
tween densities of states of the s-d exchange model 
for antiferromagnetic and ferromagnetic s-d ex- 
change coupling. This shows that incorporating 
quantum effects such as the Kondo effect is essen- 
tial for a good description of the model. 
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Method 




Features 


Simple self-consistent |9J 


0.816 


Solution of Heisenberg equation of motion, classical limit: 5 — )■ oo 


Hubbard-III [9] 


« 0.8 


Inclusion of scattering corrections 


DMFT [this work] 


^ 0.65 


Quantum calculation, local correlations 


Dual Fermions [this work[ 


^ 0.3 


Quantum calculation, fermionic non-local corrections 



Table I. Critical values of J for the metal-insulator transition in different approaches. 



The results in this work were part of the Mas- 
ter's thesis work of JS. 
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DoS's for U = 5, p=10, J positive and almost critical 



DFJ=0.35 

DMFT J=0.6 « 

DMFTJ=0.7 * 

Analytical J=0 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



Figure 1. The DoS around the critical values of J. The stars correspond to dmft results, the lines are Dual 
Fermion calculations. We see the metal-insulator taking place as the local maximum (peak) at the Fermi level 
{E = 0) disappears for higher J-values. The values are chosen such that they are just around the transition 
for the two calculational schemes, dmft and Dual Fermions. For dmft the value Jc ~ 0.65, and for the Dual 
Fermion calculation it gets as low as Jc « 0.3. 
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DoS'sforU = 5, (3=10, |J|=0.4 
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Figure 2. DoS's for |J| = 0.4. There is a difference for positive and negative J. For the local theory, dmft 
(stars), the difference is aheady visible. Beyond this approximation, incorporating non-local correlations using 
the Dual Fermion approach (lines), it is more pronounced. 



DoS's for U = 5, (3=10, |J|=0.7 



DFJ=0.7 

DFJ=-0.7 

DMFTJ=0.7 * 
DMFT J=-0.7 
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Figure 3. DoS's for \J\ = 0.7. In this regime the difference between positive and negative J is clearly visible in 
both the local (dmft, stars) and the non-local (Dual Fermion, lines) approaches. 
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